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Abstract. We evaluate the particle current flowing in steady state through a 
Bose-Einstein condensate subject to a constant force in a quasi-onedimensional 
lattice and to attractive interactions from fermionic atoms that are localized in 
various configurations inside the lattice wells. The system is treated within a Bose- 
Hubbard tight binding model by an out-of-equilibrium Green's function approach. 
A new band gap opens up when the lattice period is doubled by locating the 
fermions in alternate wells and yields an interference pattern in the transmittivity 
on varying the intensity of the driving force. The positions of the transmittivity 
minima are determined by matching the period of Bloch oscillations and the time 
for tunnelling across the band gap. Massive disorder in the distribution of the 
fermions will wash out the interference pattern, but the same period doubling 
of the lattice can be experimentally realized in a four-beam set-up. We report 
illustrative numerical results for a mixture of *^Rb and '^''K atoms in an optical 
lattice created by laser beams with a wavelength of 763 nm. 
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1. Introduction 

A quasi-onedimensional (ID) optical lattice created by the interference of two 
counterpropagating laser beams and superposed on a highly elongated magnetic trap 
provides an almost ideal periodic potential for an atomic gas. A mixture of a Bose- 
Einstein condensate and a degenerate fermion gas is an interesting system to study in 
such a confinement. The interactions can result in the spatial phase separation of the 
two components pP or in the collapse of the fermionic component [21 and generally 
induce strong changes in the atomic dynamics ||3j, depending on the values of the 
boson-boson and boson-fermion coupling constants. 

A number of experiments have examined properties of ultracold bosonic 
atoms and condensates in optical lattices: Bloch oscillations 0] Landau-Zener 
tunnelhng 0, squeezed states jZj, Josephson-like oscillations |H], a superfluid-to- 
Mott insulator transition , and ID band structure ^01 have all been observed. 
Theoretically, these effects have been studied by a variety of approaches such as the 
Gross-Pitaevski equation ^J, the discretised non- linear Schrodinger equation ^21) the 
Bogoliubov-de Gennes equations U2| i the direct diagonalization of the Bose- Hubbard 
Hamiltonian (14) . and quantum MonteCarlo simulations (15) . 

In this work we study the transport properties of the bosons in a boson-fcrmion 
mixture at zero temperature under the action of a constant driving force inside a 
ID lattice. A constant acceleration can be imparted to the boson condensate by 
using a frequency-chirped optical lattice ^7|. For a suitable choice of the system 
parameters the tunnelling of the fermionic atoms across the lattice can be neglected 
and their presence simply modifies the well depths seen by the bosons. The interesting 
cases are those in which the fermionic component introduces additional periodicities 
or generates disorder in the lattice potential. All these modifications in the lattice 
potential could be more easily realized in the laboratory by using extra laser beams 
and/or speckle ^B]. However, we choose to focus on more sofisticated experimental 
situations where the additional periodicity or the disorder is induced by the fermions, 
with the aim to study the magnitude of these effects in such a system and to address 
our numerical results to experimental groups who investigate quantum transport in 
boson-fermion mixtures. 

We describe the system by means of a ID Bose-Hubbard Hamiltonian, using a 
renormalization of the scattering lengths to embody the transverse confinement |18| . 
The Hamiltonian is not explicitly solved, but the whole lattice is reduced to a single 
dimer by a renormalization/decimation procedure Incoming and outgoing leads 
are connected to the effective dimer and the steady-state transport coefficient of the 
bosons is inferred from the evaluation of the scattered wave function of the leads due 
to the presence of the dimer. 

The paper is organized as follows. In Sec. [3 we present the physical system of 
specific interest and the model that wc adopt to evaluate on-site and hopping energies 
of the ID Bose-Hubbard Hamiltonian. In Sec. |21 we introduce a scattering matrix 
formalism for the case of out-of-equilibrium leads, with the aim of calculating the 
bosonic transmittivity through the system. The particle current is obtained from an 
extrapolation using a Landauer-like formula. The method is applied for two choices of 
the number of fermions, in cases when they are either regularly or randomly distributed 
inside the lattice. The numerical results and their physical interpretation are discussed 
in Sec.^ and the main conclusions are given in Sec.|Sl Two Appendices give technical 
details on the evaluation of some parameters entering the effective ID Hamiltonian 
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and on the recursive algorithm for the calculation of the scattering matrix, while a 
third Appendix discusses how period doubling for a pure Bose condensate can be 
achieved in a four-beam experimental set-up. 

2. The model 

We focus on a mixture of *^Rb and ''"K atoms, which has the peculiarity of having a 
large and negative boson-fermion scattering length. Let us consider the case in which 
the number Nf of (spin-polarized) fermions is much less than the number Nb of bosons 
and at most equal to the number Us of available lattice sites. With a suitable choice 
of the laser frequency ujl, the lattice potential can become considerably deeper for the 

atoms than for the ^^Rb atoms, so that the fermions can be taken as localized 
inside the lattice potential wells. 

The lattice potential for the two atomic species can be written as 

UkM^) = Rb sin'(fc^) ^ ~-^sm^{kz) (1) 

4dK,Rb 

where ^K,Rb = i^L —i^K,Rh is the detuning of the laser from the two atomic frequencies 
uk and WRb, k = 27r/A is the wave number of the laser light and r^K.Rb = dK,KbEo/h 
is the Rabi frequency, which is determined by the laser electric field Eq and by the 
atomic dipole dx.Rb- The period of the optical lattice is A/2. The atomic dipole 
depends on the natural width Fx.Rb of the transition |^ so that the ratio U^/U^^^^ is 
given by 

For the D2 lines of '"^K and ^^Rb, which are used in current experiments the 
natural widths are almost equal (Fk — F^b — 6 MHz) , while the atomic frequencies 
have values luk = 390.83 THz and = 384.227 THz. Choosing for example a laser 
wavelength A — 763 nm (or ojl = 392.88 THz) the ratio between the depths of the 
two lattice potentials is U^/U^^^ = 4.26, so that the transport of the fermions can be 
neglected in a regime where tunnelling motions of the bosons are still allowed. 

However, the presence of the localized fermions affects the transmittivity of the 
bosons by modifying the effective potential that the latter see. In the case of a negative 
mutual scattering length a well containing a fermion becomes deeper for the bosons 
inside it, so that these bosons will have a lower probability of hopping to the next 
well. 

To study the effect of the fermions on the bosonic transport as indicated above, 
we build a ID tight-binding Hamiltonian for the bosons and use a Green's function 
approach to evaluate their transport coefficient through the lattice. The Bose- Hubbard 
Hamiltonian for the bosons is 

ris 

HJ = J2[E^\^){^\+M\^){^ + l\ + \^+l){^\)] ■ (3) 

2 = 1 

Here, the parameters Ei and ji depend on the number of bosons in the lattice well 
labelled by the index i and represent site energies and hopping energies, respectively. 

We proceed to a ID reduction of the Hamiltonian by introducing the transverse 
widths (Tj_ Rb,K of the bosonic and fermionic wave functions in a cigar-shaped harmonic 
trap, together with a ID condensate wave function (l>i{z) in the i — th cell and the 
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ID fermion density jif{z). In a tight- binding scheme (t>i{z) is a Wannier function for 
the bosons in the potential U-iib{z) given in Eq.QJ and, according to the early work 
of Slater j23j . can be written as a Gaussian function, 

M^) = MO) exph(z - z.)V(2^'Rb)], (4) 
where |0i(O)p is the number of bosons in the lattice well i. The fermion density is 
similarly written as a superposition of Gaussian functions localized in a set of lattice 
wells labelled by an index i', 



nfiz) cx Vexp[(z-z,0V'^2 



(5) 



The determination of the widths o'_LRb,K and cFzRh.K is carried out variationally, as 
described in [Appendix A| 

We can now evaluate the parameters entering the effective Hamiltonian. The site 
energies are given by 



/i^v^ 1 
-7; 1- Uiih{z) + 7:5b&l</'i(^)P + 9bfnf{z) - mRhaz + Cnb 



(6) 



where mRb is the *^Rb mass, a = F /niYih is the acceleration due to a constant external 
force F acting on the bosons, Cnh — / {2mnh(y\Yih) + ^"^Rb^iRb'^'j.Rb' '^J-Rb 
is the radial frequency of the harmonic trapping potential. The parameters gbb and 
Qhf are the strengths of the ID boson-boson and boson-fermion interactions, which 
are given by 

- 47rfe^ abb 
9bf = - 



I 1 , 1 \ 

with abb, ctbf the boson-boson and boson-fermion scattering lengths and rUr the 
boson-fermion reduced mass. Consistently with the model, the hopping energies 7^ 
are given by 

Unhiz) + -gbb\(t)iiz)f + gbfUfiz) + Cnb 



7i = / dz(t)i{z) 



Mi{z)- (8) 



2TORb 2- 

This completes the determination of the effective ID Hamiltonian for the bosons, that 
we shall use to evaluate the transmittivity for various distributions of the fermions in 
the lattice wells. 



3. The transmittivity 




Figure 1. Schematic representation of the tight-binding Hamiltonian for a 
condensate in a ID lattice, with some sites containing a fermion. 

In Fig. n we show a schematic representation of the Hamiltonian Q in the absence of 
an external force. For such a system the boson site energies can be of two types, Ei 
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and E2 depending on the absence or presence of a fermion in the weU, and the boson 
hopping energies can take three values: 71 if the hopping is between a site with a 
fermion and a site without, 72 in the case of two neighboring wells which both contain 
a fermion, and 73 in the case of two neighboring wells without fermions. We shall 
consider the following situations: (i) the fermionic band is empty (EF), namely no 
fermions are present in the lattice; (ii) the fermionic band is filled (FF), namely there 
is one fermion per well; (iii) the fermionic band is half-filled {Nj = ns/2) with every 
other well occupied in an orderly way (OHF); and (iv) the fermionic band is half-filled 
but the fermions are distributed on the lattice in a disordered way (DHF). 

In the tight-binding regime the system of bosons behaves in the first two cases 
as a one-band system, which is described by a single site energy {Ei in the first case 
and E2 in the second) and by a single hopping energy (73 and 72, respectively). The 
third case describes a two-band boson system with two alternating site energies (Ei 
and E2) and a single hopping energy (71). There is in this case a folding of the 
dispersion relation in quasi-momentum space, with the opening of an energy minigap 
AE = \E2 — Ei\ between two sub-bands of equal energy width. In the fourth case the 
bosons again form a two-band system of the type illustrated in Fig. ^ and the presence 
of disorder will cause a broadening of the two bands. We show in [Appendix B| how 
the third case could be realized in the laboratory by using a four-beam set-up instead 
of an ordered array of fermionic atoms at half occupancy. The forth case could be 
realized by adding some speckle to the lattice |16j . 




Nout 



Figure 2. Band tilting due to a constant force F and coupling of the condensate 
to an incoming and an outgoing lead for the case of (a) a one-band system and 
(b) a two-band system. The position z is in units of A/2. 



The effect of applying a constant force to the condensate can be described as a 
tilting of the bands in coordinate space, due to the induced position dependence of the 
site energies. This is shown in Fig. 12 the situation is analogous to that obtained in 
an electron gas on a lattice by applying to it a static, uniform electric field [SI]- The 
transmittivity of such a system can be calculated, following the proposal of Farchioni 
et al. |25j , through an evaluation of the scattering matrix which explicitly takes into 
account the presence of the bias [221 ■ ^'^^ this purpose, let us connect the system to 
incoming and outgoing leads which mimic its coupling to the continuum by injecting 
and extracting a steady-state particle current, respectively. The bosons are initially 
in their lowest momentum state {p = 0) and can leave the lattice only at p = ft,/ A (in 
the one-band case) or at p = (in a small-minigap two-band system) if the coupling 
to the leads is effected as is shown in Fig. |21 Elastic Bragg scattering can also occur 
at both these points in quasi-momentum space as well as at the top of the lower sub- 
band in the split-band configuration. As we shall see, in the latter case this produces 
resonances in the transmittivity. 
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Let us define iVout (with iVout ^ '^s) a-s the number of hops after which the 
condensate reaches the highest energy point in the dispersion curve and can either 
leave the lattice or be Bragg reflected. The two leads have a band width and their 
centers at E™ and E'™* are shifted by F7VoutA/2 in order to optimize the couphng 
with the lattice. The outgoing wavefunction (j)^^^{z) is an eigenstate of the outgoing 
lead and is connected to the incoming wavefunction (j)'^^{z), which is an eigenstate of 
the incoming lead, by the relation 

C,t(z) = f^+(£;)CW- (9) 
The wavevectors k and k of the two wavcfunctions arc uniquely identified by the 
relation 

~ 2 

K,K = — arccos . 10 

X \ 2jL J 

The operator n+{E) = I + G°{E)T{E) in Eq. © involves the Green's function 
G^{E) = {E — Hq)~^ of the two leads, as well as the scattering matrix T{E) ~ 
Hi\i — G^{E)Hi\~^ where Hj is the Hamiltonian acting on the first A'out sites of the 
lattice subject to the force F (see [Appendix C| |. 

In the next step we renormalizc the interaction Hamiltonian Hi into an effective 
Hamiltonian Hj{E) acting in the subspace {| 1) , | A*"out)} formed by the two edge 
sites that are directly connected to the leads. In this renormalization approach the 
scattering matrix is also projected on the subspace {| 1) , | -/Vout)} and can be written 
as a 2 X 2 matrix and evaluated with a recursive algorithm as described in |Appendix C| 
In such a formalism the transmittivity coefficient T{E), which is defined as 

^/^N ^ hm„^+oc {n\ Cm ) (Cut I ") ^out ,-|^^. 
lim™^_oo(m|Cn)(CJnT-)win 
with Win and Wout the velocities of the incoming and outgoing wavefunction, is related 
to the off-diagonal element T{E)i^N^^t of the scattering matrix by 

sm (^K-jsm [n-^ 

It is easy to show that Eq. (|12|l is equivalent to the trace of Eq. (63) in the work of 
Paulsson 26 . 

The transmittivity determines the particle current passing through the lattice, 
and to calculate it we can exploit Landauer's approach adapted to a bosonic gas with 
out-of-equilibrium leads. Again following Paulsson |26| . the current j can be written 

J = ^J - Ml) - fiE mWiE) (13) 

where f{E— fi) is the Bose-Einstein distribution and /^i,2 are the chemical potentials 
of the reservoirs which act as the source and the sink of the bosons. Note that the 
normalization Nh/us of the bosonic wavefunction plays the role of the electric charge 
in the Landauer formula |27| . In the zero-temperature limit we have 

f{E Ml) - f{E m) ^ ^^^Qj^ (M2 - A^i) (14) 

and the density probability df{E — jj)/ BE is a delta function 5{E — fj). Then Eq. ({T^ 
gives 

J = -^T{jl)Af, (15) 



r(i^)^4 '-^- sin (.9) sin (^9). (12) 
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where A/i = M2 ^ /^i = is the difTerence in chemical potentials and 

JI = (/ii + /i2)/2 is their average value, corresponding to the dashed line in Fig. [3 

The transmittivity may also be expressed in terms of the number of bosons exiting 
from the lattice as a drop emitted under a constant driving force, as in the experiment 
of Anderson and Kasevich 5 . The current is written in that context as 

j = Ndropl'B, (16) 

where fB = {F^)/{'^h) is the frequency of Bloch oscillations and iVdrop is the number 
of bosons in the first drop. This can be calculated from Eas. (|15ll6ll as 

%^ = ^n-p). (17) 

Equation lfT7|l holds for a steady-state current and will not be applicable to the 
subsequent drops unless the condensate is continuously replenished. 



4. Results 

We have calculated the transmittivity of the model as a function of the acceleration a 
for the four cases listed at the beginning of Scc.O The difference A/i of the chemical 
potentials of the two leads does not depend on a and is fixed by the total band width 
47L of t he le ads, so that the current j is directly proportional to the transmittivity 
(see Eq. (|15|) '). We have chosen to focus on the typical experimental parameters used 
at LENS for the ^^Rb-^^K mixture, setting TVf, = 10^ n, = 200, U^^ = 3.5£;rRb, 
and w^Rb = 27r X 85.7 Hz. As previously remarked, we set A = 763 nm so that we can 
consider each atom as localized in a lattice well. Moreover, for this choice of the 
laser wavelength, the detunings (5k, Rb are much larger than the natural widths rK,Rb 
and spontaneous emission can be neglected. 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 



Figure 3. The transmittivity as a function of the acceleration a (in units of the 
acceleration of gravity g) for the cases of an empty fermion band (EF: +) and a 
completely filled fermion band (FF: X). The lines are guides to the eyes. 

In Fig.0we show the transmittivity of a one-band system as a function of a for the 
one-band cases in which either the fermions are absent (EF) or each well is occupied 
by one fermion (FF). The effect of the presence of the fermions should be that the 
boson- fermion attractions diminish the hopping probability of the bosons between two 
neighboring wells, so that the transmittivity should be lowered. However, for the case 
of only one fermion per well and for the usual values of the boson-fermion scattering 
length, this effect is not visible (see Fig. O . 
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Figure 4. (a) The transmittivity at lialf occupancy as a function of the 
acceleration a (in units of the acceleration of gravity g) for an ideal two-band 
system (OHF: +) and for randomly distributed fermions (DHF: X). (b) The same 
transmittivity data are plotted against the ratio Tg/r. The lines are guides to 
the eyes. 



In Fig.01we show the transmittivity of the two-band system for Nf — 100 fermions 
which are distributed in the lattice wells either regularly (OHF) or completely at 
random (DHF). The value of A^out in these calculations has been taken to run from 
200 to 12: this corresponds to going from left to right in Fig.0fa) and from right to 
left in Fig.0fb). The transmittivity in Fig. 01 is on average strongly reduced relative 
to the cases shown in Fig. |3 but in the OHF case it shows a very prominent structure 
of peaks and troughs. The observed pattern arises from the interference occuring at 
p = -~h/2X (the left-hand edge of the Brillouin zone in the split-band configuration) 
between the boson wave packets that are Bragg scattered a.t p — h/2X in the lower 
sub-band and those that have tunnelled into the upper sub-band, travelled through 
it, and tunnelled back into the lower sub-band at p = —h/2X. The two times that 
govern the interference pattern are the Bloch oscillation period Tg = h/{FX) and 
the time t taken by the bosons in coherently tunnelling twice across the minigap, 
which is proportional to (HNb/AEns) |2H1- We expect constructive interference in 
the reflectivity at p = in the upper sub-band, i.e. destructive interference in the 
transmittivity, whenever the ratio Tb/t is an integer number. This behaviour is 
illustrated in Fig. 0fb) by replotting the transmittivity data against Tb/t with the 
choice T = {Stt'^ /8){hNi,/ AE Us). Of course this specific value of the proportionality 
constant in the tunnelling time, which is consistent with the Heisenberg principle, 
depends on the model Hamiltonian that we have assumed. 

Returning to Fig.0fa), the effect of disorder in the distribution of the fermions in 
the half-occupancy DHF case is to strongly diminish the contrast in the interference 
pattern. The disordered distribution that we have examined in Fig. 01 has been 
produced from a random number generator and the pattern is completely washed 
out in this case of maximum disorder. While the average occupancy of a site by a 
fermion is only one-half, the transmittivity of the bosons is greatly reduced relative to 
the FF case in Fig. O Evidently, the main effect at play is due to the fluctuations in 
the disordered distribution of the fermions, which tend to break the coherence of the 
Bose condensate. 
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5. Concluding remarks 

In conclusion, we have studied the condensate current induced by a constant external 
force in a ID optical lattice in the presence of various numbers of fermionic atoms, 
distributed over the lattice in orderly or random ways. We have focused on the ^^Rb- 

mixture, with a choice of the laser generating the lattice potential such that the 
barriers for the fermions are much higher than those for the Bose condensate. In this 
regime fermionic transport can be neglected, but the presence of the fermions and their 
spatial distribution affect the bosonic transport via mean-field attractive interactions. 

We have calculated the transmittivity and the condensate current without 
explicitly diagonalizing the Bose-Hubbard Hamiltonian, but by exploiting an out-of- 
equilibrium Green's function approach adapted to the bosonic case. When there is one 
fermion per well, the current is essentially the same as in the absence of fermions and 
increases monotonically with the driving force. When every second well is occupied by 
a fermion, instead, the energy spectrum breaks into two sub-bands and the probability 
that the condensate may tunnel through the band minigap depends on the energy 
that it has acquired at the edge of the first sub-band and on the energy width of the 
minigap. Part of the condensate is back-scattered when it arrives at the edge of the 
lower sub-band and can interfere at the opposite edge of the Brillouin zone with the 
part that has tunnelled to the upper sub-band and travelled through it. This gives rise 
to an interference pattern showing marked maxima and minima in the transmittivity, 
whose location is an integer multiple of a parameter directly proportional to the energy 
minigap and inversely proportional to the driving force. 

The interference pattern loses resolution when the distribution of the fermions in 
the lattice becomes disordered , up to being completely washed out in the case of full 
randomness. An ordered period doubling of the lattice could be realized, however, 
with a four-beam experimental set-up. Observation of resonances in the current as a 
function of the driving force would allow a measurement of the tunnelling time through 
the minigap. 
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Appendix A. Determination of radial and Eixial widths 

We adopt the variational method proposed for Bose condensates by Salasnich et 
al. |18j . Within a ID nonpolynomial nonlinear Schrodinger equation approach, we 
find that the evaluation of the transverse widths a±_ Rb,K requires the self-consistent 
solution of the Euler-Lagrange equations 

Assuming that ctj^ Rb,K are essentially constant along the axial position, we set 
in Eas. HA.l|) |0j(2:)P ~ Nhl{\fT^nsOz'Rb) and ~ l/(\/7ra-^K) or = 0, 



(A.l) 
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depending on the presence or absence of a fermion in the corresponding weU. The 
axial widths cr^ R,b,K entering the axial densities depend themselves on the radial widths 
f± Rb,K- Consequently we first solve Eas. (|A.ip for given values of azRh,K, evaluated in 
the absence of interactions and then, having calculated the radial widths, we evaluate 
the correction to the axial widths. 

Minimization of the bosonic and fermionic energies with respect to ctz Rb.K yields 
the two coupled equations 

. lm„u,,2 _2 I 9bbNb I 9b f Cr^Rb 

ZmnhCTzBh ^ ^"iRbCT^Rb 2V27rnscr^Rb V (cr^ + k) 

9b!Nb '2 



'^JZ—p. ^"iKt^zK'^zK + 7Z2 ,2 N3/2 - ^ 



for the case of a well containing a fermion, while in the absence of a fermion we have 

1_ 2 2 , 9bbNb 



^"T,RbCT^Rb ^ 2V27rnsCTzRb 
as already found by Chiofalo et al. |13|. Here, the axial frequencies WzRb.K for a single 
well are expressed in terms of the potential depths U^^ ^ and of the recoil energies 
-ErRb.K = ft^fc2/2TORb,K as 



t^zRb.K = l:^/4t/{^b,K^''Rb,K■ (A.4) 



The solution of Eqs. HA.l|) - (jA.3|l provides the further input that is needed in Eqs. © 
and ©. 



Appendix B. Period doubling with a four-beam set-up 

Four laser beams can be used to create a ID lattice with a double periodicity along the 
z direction. Two beams are set in a counterpropagating configuration along the z axis, 
while the other two beams are rotated by angles of 60° and 120° with respect to the 
z axis. In the region of space where the four beams are superposed, the longitudinal 
electric field E is proportional to 

E (X /3e"(e*''^/2 + g^'^ie^^'^^/^) ^ ^^^^kz ^ ^i4,^^~ikzy -^^ 

In Eq. HIj.l() (3 and a are the relative amplitude and phase of the rotated beams with 
respect to the counterpropagating pair, is the relative phase of the rotated beams, 
and (/)2 the relative phase of the counterpropagating beams. All these parameters can 
be set by choosing suitable optical components. 

Upon setting (pi = (j)2 — 'n and a — 7r/2, the resulting lattice potential will have 
the shape 

U{z)o^0^sm^{kz/2) + s\T?{kz) (B.2) 

and the double period enters through the well depths, that is the site energies, which 
will differ by a factor Instead, by setting = 02 = and a — 7r/2 we have 

U{z)<^0^cos^{kz/2) + cos^{kz) (B.3) 

and the potential hills (that is the hopping energies) have alternating values. In the 
case of a ^^Rb condensate, a lattice potential of the type in Eq. (|lj.2() with /3 ~ 0.1 is 
equivalent to the effective potential generated by a regular distribution of '^'^K atoms 
inside every other well of an optical lattice created with a single pair of laser beams. 
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Appendix C. The recursive algorithm 

Figure ICTI shows a schematic representation of the whole Hamiltonian for the system 
composed of a finite lattice and two leads, and of its decomposition into the sum of 
Hq and Hj for the evaluation of the scattering matrix. We show below how Hj can 
be reduced to an effective Hamiltonian Hj for a dimer in a subspace {| 1), | A^out)} 
and how this leads to an expressions for the scattering matrix and for the outgoing 
wavefunctions . 




ll> IN„„,> 



Figure CI. Schematic representation of the Hamiltonian H = Hq + Hj of the 
system composed by the leads (-ffo) a-^d by the lattice {Hj), and of the effective 
Ifamiltonian Hj projected on the two edge sites. 

Let n be the site index number. To calculate the transmitted wavefunction we 
have to evaluate 

{n I Cut) =Y.^n\l + G\E)T{E) \ m) {m |C) (CI) 

for n > iVout- Since the wavefunctions (j)'^^ and 0^^^ are defined in two disconnected 
spaces, Eq. (|C.1|I becomes 

(^ICut)= E {n\G^{E)\l){l\T{E)\m){m\<l^l) 

m<l,l 

= {n I G°(i?) I iVout)(iVout I TiE) 1 1)(1 IC) 

= V2[G°(S)]„,w„„, [T{E)]n^^^^i sm{K\/2)u,{X/2) , (C.2) 

where we have set (1 l^fjj) = MK(A/2)[exp (iKA/2) — exp (— iKA/2)]/iV2 with Uk(A/2) = 
4>i{z). It can be shown that Eq. (|C.2p is equivalent to Eq. (34) in the work of 
Paulsson ^ by making use of the relation G°{E)T{E) = G{E)Hi{E). 

The Green's function element [G°(£')]„,jVout Ea. (|C.2|l determines the coherence 
between site n and site A^^out on the chain for the outgoing lead. It can be written as 
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in terms of quantities defined in the main text. The scattering matrix element 
[T{E)]]y^^^ i can be evaluated from the effective Hamiltonian Hi{E) as schematically 
illustrated in Fig. ICll 

[TmN...,i = {Nout \HiiE)[h - &iE)HiiE)]-'\l) (C.4) 

where I2 is the 2x2 identity matrix and G^{E) is the Green's function of the 
leads projected onto the subspace {| 1), | Na^it)}- The only non-zero elements of the 
2x2 matrix &{E) are the diagonal ones, which are given by Q{E — E™,^]^) and 
g{E-El^\-fL), with 

1 



g{E -EsnL) = ^[E-Es-J{E-EsY- AjI] . (C.5) 

Finally, the elements of the matrix Hi{E) in Eq. (jC.4p can be calculated by 
decimation of the sites (2,3,..., A'out — 1) and can be written in the form of continued 
fractions. That is, 

[Hj{E)],^, = [Hj],s + (C.6) 

E - [i?/]2.2 - TTTy^ 

and 



1 , 1 



[^^/(^^)]i,^o„. - [g/]i,2 ^ [gj]2,3 [Hih.A . . . , (C.7) 

E - [tll\2,2 ^ ^ Wl\^^2 

E — [-ff/]2,2 



with [Hi{E)\n,^^,i = [Hi{E)]i^M„^i and [Hi{E)]N,^^,N^^i given similarly to [Hi{E)]i^x 
in Eq. iTTfill . 
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